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Abstract 

We establish local balance equations for smooth functions of the vorticity in the DiPerna- 
Majda weak solutions of 2D incompressible Euler, analogous to the balance proved by 
Duchon and Robert for kinetic energy in 3D. The anomalous term or defect distribution 
therein corresponds to the "enstrophy cascade" of 2D turbulence. It is used to define 
a rather natural notion of "dissipative Euler solution" in 2D. However, we show that the 
DiPcrna-Majda solutions with vorticity in LP for p > 2 are conservative and have zero defect. 
Instead, we must seek an alternative approach to dissipative solutions in 2D. If we assume 
an upper bound on the energy spectrum of 2D incompressible Navier-Stokes solutions by 
the Kraichnan-Batchclor k~ 3 spectrum, uniformly for high Reynolds number, then we show 
that the zero viscosity limits of the Navier-Stokes solutions exist, with vorticities in the 
zero-index Besov space B^'°° ', and that these give a weak solution of the 2D incompressible 
Euler equations. We conjecture that for this class of weak solutions enstrophy dissipation 
may indeed occur, in a sense which is made precise. 
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1 Introduction 

In 2-dimensional turbulence it is the enstrophy Q(i) := ^||cc(t)||2 that is expected to cascade to 
small length-scales, not the energy as in three space dimensions |lJ-(3[]- In a view that goes back 
to Onsager jij , such turbulent cascades are conjectured to be described, in the limit of infinite 
Reynolds number, by singular (or weak) solutions of the incompressible Euler equations. More 
recently, Duchon and Robert H have shown how Onsager's idea of a dissipative Euler solution 
may be formalized in the three-dimensional case via a local energy balance relation. It is our 
purpose here to similarly formalize the notion of a 2-dimensional dissipative Euler solution, 
corresponding to the enstrophy cascade. 

We consider weak solutions of the 2D Euler equations in the vorticity- velocity formulation: 

d t u + (u-V)w = 0, (1.1) 

with u = K *uj given by the Biot-Savart kernel K. We show first that when the vorticity fields 
u;(x, t) are suitable measurable functions and (1.1) is interpreted in the sense of distributions, 



then a local balance is satisfied 

d t h{uj) + V-[uh(u)} = -Z h (co), (1.2) 
for nonnegative, convex functions h(uS). Of course, (1.1) formally just expresses the conservation 



of vorticity along fluid particle trajectories, so that ( |1.2| ) would naively be expected to hold with 
Zh(uj) = 0. Dissipative weak solutions might be taken to be those for which this distribution 
is nonnegative: Z^uj) > 0. The balance equation ( |1.2j ) makes more precise Polyakov's analogy 
of the enstrophy cascade with conservation law anomalies in quantum field-theory (such as the 
axial anomaly in QED) ||. The distribution appearing as a sink term on the right side of ( |1.2[ ) 
corresponds closely to such an anomaly. However, we show under rather general conditions, 
even weaker than those in our earlier work |?J, that Zh(u) = 0. For example, we show that the 
anomaly vanishes for functions h of "power-p growth" at large arguments, whenever the initial 
vorticity satisfies an L p bound in space. In particular, this means that enstrophy is conserved 



by a 2D Euler solution whenever the enstrophy itself is finite. [] This state of affairs presents a 
striking contrast with the situation in 3D where, as discussed by Duchon and Robert ||], energy 
dissipation is expected to be possible for incompressible Euler solutions with finite energy. 

The above results necessitate an approach to the notion of dissipative Euler solution in the 
2D case which is therefore rather different from that of Duchon-Robert for 3D. Nevertheless, 2D 
turbulence theory is still a useful guide to the correct formulation. Indeed, the above results are 
in perfect agreement with the expectations of the classical theories which predict that 

the small-scale energy spectrum in the enstrophy cascade range of 2D turbulence shall be of 
the form E(k) ~ Ck~ 3 (with at most a logarithmic correction). Hence, the classical theories of 
2D turbulence predict an infinite total enstrophy but a finite spectral flux of enstrophy. In the 
following we shall formalize this notion for an appropriate class of weak Euler solutions in 2D. 
We define as dissipative those solutions which have a nonnegative flux of enstrophy (possibly 
zero or infinite) asymptotically to infinitely high wavenumbers. The relevant solutions must, 
however, have vorticity fields which exist only as distributions and not as ordinary (measurable) 
functions. We show that such solutions of 2D Euler equations with a Kraichnan-Batchelor A; -3 
energy spectrum are obtained as zero-viscosity limits of the Leray solutions of 2D Navier- 
Stokes, whenever upper bounds on the spectrum by the Kraichnan-Batchelor prediction hold 
uniformly in the viscosity. We then show that the notion of enstrophy flux is well-defined 
for such distributional solutions of 2D Euler, although the enstrophy itself may diverge. Our 
natural conjecture is that the flux is asymptotically nonnegative at small length-scales for all 
such "viscosity solutions" and, for suitable initial data, even strictly positive. 

Our main results are stated as three Theorems in the following Section 2, where their content 
is further discussed in detail. The proofs of the Theorems are outlined in the final Section 3. 

This statement has, among other consequences, the implication that no power-law 2D energy spectrum 
E(k) ~ Ck~ n with n > 3 may be compatible with existence of an enstrophy cascade. The conformal "solutions" 
studied by Polyakov Q that have spectral exponent n > 3 therefore cannot exhibit an anomaly in the enstrophy 
conservation law, as he has proposed. 
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2 Statement of Results 



Before stating precisely our theorems, it will help to motivate the statements (and the proofs) to 
give a brief, heuristic argument for the existence of the enstrophy cascade. In |?J we considered 
a "filtered" form of the 2-D Euler equations (see also ||): 

dtu £ + V-[u £ oj £ + a £ ] = 0, (2.1) 

where ui £ = (p £ * u> for a smooth mollifier ip, <£> e (x) = e~ 2 </?(e x), and <y £ = (ucj) e — u £ uj £ . 
The new term a £ represents a turbulent spatial transport of vorticity due to the eliminated 
small-scales. It is straightforward to show that the balance holds that 

d t h(uj £ ) + V-[u £ h(uj £ ) + ti{u £ )cT £ ] = ti'(u e )Vuj e -<T £ . (2.2) 

The term Zh^ £ (uj) := —h"(uj £ )'Vuj £ -(T £ represents a transfer of /i-stuff from length-scales > e to 
smaller scales. Based upon the notion of "UV-locality of interactions", a natural approximation 
is to take cr £ « (const.) [(u £ lo £ ) £ — u. £ lj £ } and then to Taylor expand to leading non- vanishing 
order to obtain 

a £ « Cs 2 B £ -Voj £ . (2.3) 

Here D £ is the filtered velocity-gradient tensor Dij = dui/dxj; also, a spherically symmetric 
mollifier has been assumed. The first of our approximations is analogous to the "similarity 
model" employed by engineers in large-eddy simulation of three-dimensional turbulence and the 
second to its further simplification, the "nonlinear model" ||. The matrix D £ is traceless and 
has, in vortical regions of the flow, a pair of imaginary eigenvalues and, in strain-dominated 
regions, two real eigenvalues of equal magnitude S £ but opposite signs. It stands to reason 
that, in the latter straining regions, the compression of vorticity level sets will tend to align 
the direction of the vorticity gradient Vu; £ with the eigendirection of D e corresponding to the 
negative eigenvalue. Indeed, such alignment has been observed in simulations to hold (for the 
unfiltered quantities) with a high probability Assuming it to hold exactly, we find that 

a £ « -Ce 2 S £ Vu £ . (2.4) 
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This is precisely an eddy-viscosity model, with effective viscosity v £ = Ce 2 S £ at scale e. It leads 
to an effective dissipation Z £ (uS) ~ — v £ \^u> £ \ 2 . If the vorticity field is Holder continuous with 
exponent s, uj G C s , then Vu> £ ~ e s_1 for small e and S £ ~ S independent of e. In that case, 
Z £ {oS) ~ e 2s for e — > 0, so that we expect an asymptotic enstrophy cascade only when s = 0. 
This is precisely the "mean-field" scaling exponent in the Batchelor-Kraichnan theory jjj]-|H]- 
We now state our main theorems: 

Our first theorem establishes the local vorticity balance equations for the weak Euler solu- 
tions constructed by DiPerna and Majda for initial data luq £ L p ,p > 1 [O]. Although they 
considered solutions in the whole plane M 2 , we shall restrict attention for simplicity to solu- 
tions on the 2-D torus T 2 . DiPerna and Majda also established existence of weak solutions in 
the velocity-pressure formulation, but it is not hard to show that, for p > 4/3, the associated 
vorticity field in their solution also satisfies the weak vorticity- velocity equations (see below). 
In fact, the only property of the DiPerna-Majda solution that we will employ in our proof is 
that to G L°°([0, T], L P (T 2 )) and our theorem would apply to any other such solutions as well. 
To state our theorem, we must introduce an appropriate class of differentiable functions 

H p := {h\ h e C^R), \h'(u)\ < C(l + M^ 1 ) for some C > 0} (2.5) 

which have at most L p -growth. We then have the following: 

Theorem 1 If u G L°°([0, T], L P (T 2 )) and the associated u = K * u for p > 4/3 are a weak 
solution of 2-D incompressible Euler in the vorticity-velocity formulation, then for h G Tt r f] C 2 , 
with r = t;P — 1 for | < p < 2, r < p for p = 2, and r = p for p > 2, the balance holds 

d t h(Lo)+V-[uh(uj)} = -Z h (oj) 

in the sense of distributions. The righthand side is given by the distributional limit 

Z h {u) = lim -h"{uj £ )Vuj £ -a e (2.6) 

which exists for any choice of mollifier ip which is C°° , nonnegative, and compactly supported, 
with unit integral, and it is independent of that choice. For the special case of the enstrophy 
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integral, h(u>) = ^\u>\ 2 , when p > 2, we write simply Z(u>) = Z^iuj). In that case, there is the 
alternative expression: 

Z{uj) = lim- f d 2 £ V^ £ (£)-A^u|A^| 2 (2.7) 
e-^o 4 J 

where A^o;(x, t) = u(x + £, t) — uj(x, t), likewise for A^u, and ip is further restricted to be an 
even function of its argument. 

Note that, formally, Zh(uj) = h"(oj)Z(uj), so the fluxes of general convex functions are, in some 
sense, proportional to the enstrophy flux with a nonnegative factor. The last expression ( |2,7D 
for the enstrophy flux has a nice interpretation as a local, non-ensemble-averaged form of the 
"-2 law" for the direct cascade, in its form applicable without isotropy (see @, Appendix B). 
Thus, the defect distribution in the vorticity balance equations has an exact connection with 
the enstrophy cascade in 2D turbulence theory. 

However, we next show that this distribution is, in fact, zero for the DiPerna-Majda weak 
solutions, which therefore conserve the integral 

I h (t) = J d 2 x h(u(x,t)) (2.8) 

for all h of suitable growth: 

Theorem 2 If to G L°°(0, T; L P (T 2 )) is a DiPerna-Majda weak Euler solution for p > 2, then 

d t h(u) + V-[uh(u)] = (2.9) 
in distribution sense for all h € 7i r with r = p when p > 2 and for any r < p when p = 2. 

In Q it was proved that such a conservation statement holds for uj € L p (0, T; 5p'°°(T 2 )) for 
s > 0,p > 3 where Bp°°(T 2 ) is the standard Besov space of functions in L P (T 2 ) which are Holder 
of index s in the L p -mean sense That theorem was thus analogous to the Besov-space 

improvement of Onsager's original conservation result for 3D, which was proved by Constantin, 
E, and Titi jOJ. We now see that the smoothness assumed in |7| was unnecessary and that 
simple LP bounds alone are sufficient for conservation. Essentially the same result was already 
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obtained by P.-L. Lions in [13], Section 4.1, based upon his earlier work with R. J. DiPerna |15[ . 
He showed there that the DiPerna-Majda solutions with p > 2 are "renormalized solutions" in 
the sense of DiPerna-Lions [15|, which amounts to the requirement that ( |2,9D hold. In fact, 
global conservation 

/ d 2 x h(u(x,t)) = I (f 2 x /i(£j (x)), t>0 (2.10) 

is shown in |l4| to hold for all h £ W p even when p = 2, just as in the proof of Theorem II. 2 
and equation (26) in DiPerna-Lions jy|. In particular, taking h{uj) = ||w| 2 , a remarkable 
statement is true that enstrophy dissipation is not possible for any 2D Euler solutions with finite 
enstrophy. We conclude more generally that the DiPerna-Majda weak solutions are not relevant 
to the problem of constructing dissipative Euler solutions. In the language of turbulence theory, 
they do not support enstrophy cascades over infinitely-long ranges of wavenumber. 

The conservation properties of the DiPerna-Majda solutions for p > 2 have an intuitive 
explanation. It has been noted recently that breakdown of uniqueness of Lagrangian particle 
trajectories in Holder but non-Lipschitz flows can be a mechanism for the anomalous dissi- 
pation of the analogous integrals as ( |2.8| ) for passive scalars Jl7]]-|l!J. For the 3D problem, 
Shnirelman has found a weak solution which dissipates energy globally, by constructing a gen- 
eralized flow with random Lagrangian trajectories [20]. In the case of the Yudovich solutions 
of 2D Euler with u £ L°°(T 2 ) [21], it has long been known that they are conservative pre- 
cisely because the corresponding velocity field is log-Lipschitz and the Lagrangian flow maps 
Xt are unique, volume-preserving homeomorphisms. Therefore, the Yudovich solution is given 
simply by cj(x,t) = u;o(^-i( x )) i n terms of the inverse-Lagrangian map. All of the integrals 
Ih(t) m ( |2.g| ) are then trivially time-invariant. DiPerna and Lions in their paper |15] show 
that there are likewise unique Lagrangian flow maps Xt(x) with X € C(0, T; L P (T 2 )) whenever 
u € -^ 1 (0, T; W 1,P (T 2 )) for p > 1 and that these maps preserve Lebesgue measure when V-u = 0. 
The "renormalized" solutions of the linear advection equation constructed by DiPerna-Lions 
2 The same remark was made in a recent preprint of E and Vanden-Eijnden 
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are shown to have precisely the form w(x, t) = uJo(X^t( x ))- While it is not true in general that 
the distributional solutions of 2D Euler in the sense considered here are renormalized solutions, 
Theorem 2, as we have observed above, shows that this is so for the DiPerna-Majda solutions 
when p > 2. Hence, the conservation properties of these solutions are again connected with the 
uniqueness of Lagrangian particle trajectories. 

The above results are negative — in the sense that they imply a lack of enstrophy dissipation — 
but we wish to emphasize that they are fully consistent with the expectations of 2D turbulence 
theory. In fact, the Navier-Stokes solutions exhibiting an enstrophy cascade are expected to 
have the Batchelor-Kraichnan energy spectrum 

E(k,t) ~ C?7 2/3 (t)A;- 3 , (2.11) 

where r](t) is the enstrophy dissipation rate per volume (l|, This spectrum should hold 
at high wavenumbers k 3> ko(t), the wavenumber of peak enstrophy, up to a wavenumber 
kd{t) = i^~ 1//2 77 1 / 6 (i), at which the dissipation by viscosity v becomes relevant. Equivalently, 
the enstrophy spectrum predicted by Batchelor-Kraichnan theory is 

fi(M) ~ C77 2/3 (i)Ar\ (2.12) 

for ko(t) <C k <C kd(t). In the limit as v — ► this spectrum extends all the way to +oo and 
its integral diverges, implying an infinite total enstrophy. As we have seen, this is rigorously 
required to have limiting Euler solutions which can dissipate enstrophy. 

In fact, velocity fields u(t) with the Batchelor-Kraichnan spectrum ( [2.11D for all k 3> ko(t), 
when that spectrum is interpreted in a suitable sense, must consist of u(t) £ B^ 00 ^ 2 ), with 
corresponding vorticity u(t) € B®' 00 ^ 2 ), the Besov space of zero index. The definition of spec- 
trum which is relevant is a "Littlewood-Paley spectrum" which was earlier used by P. Constantin 
in [^U to prove a rigorous upper bound. This spectrum is defined in terms of the Littlewood- 
Paley decomposition of the velocity u(t) = ]C!/v=o u N(t) with ujy(t) = -0jy * u(i), for a smooth 
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partition of unity in wavenumber space 

oo 

Mk) + ^ $ N (k) = 1 (2.13) 

N=l 

where supp(V^v) C [2 N ~ 1 ,2 N+1 ] for N > 1 and supp(^ ) C [0,2]. The Littlewood-Paley 
spectrum is then defined by 

E LP (k,t) := k^WuN^fv (2.14) 

for k G [2^, 2 Ar+1 ). With this definition it is not hard to see that u(i) £ £? 2 ,QO (T 2 ) precisely when 
the spectrum satisfies a bound of the form Ej J p(k,t) = 0(k~ 3 ). In fact, the Littlewood-Paley 
criterion for / € -Bp'°°(T 2 ) with — oo < s < oo and p > is just that 

j.,00 : = sup 2^11/^1 1 lp (2.15) 

P N>0 



be finite, and for p > 1 this is a norm making .Bp' (T ) into a Banach space. See [12|, Sections 



2.3.1-3. Note we assume only big-0 bounds on the spectrum and not a power-law scaling. In 
fact, not long after his first paper on 2D turbulence, Kraichnan argued that there should be 
a logarithmic correction, E(k,t) ~ Cr] 2 / 3 (t)k~ 3 [\n(k/ko)}~3 , where ko is the lower end of the 
enstrophy cascade range [|23[|. In any case, it is still generally believed that the true energy 
spectrum must be bounded above by the form (2.11) at high Reynolds number. In that case, 



we see that the vorticity field u)(t) has the Besov index s = but not necessarily any larger 
index. As we have already remarked, u(t) £ Bp°°(T 2 ) with p > 2 and s > could not be 
consistent with a non-vanishing enstrophy dissipation. 

It is still an open question whether solutions of 2D Euler equations exist with velocities and 
vorticities in such Besov spaces and, if so, whether they dissipate enstrophy in a suitable sense. 
We shall prove here a few simple results in this direction and, in particular, advance our main 
conjecture. 

We show first that an upper bound on the energy spectrum of the solutions of the 2D 
Navier-Stokes solutions u u (t) by the Batchelor-Kraichnan spectrum, when that bound is uni- 
form in the viscosity, implies the existence of 2D Euler solutions u(i) in the appropriate Besov 



spaces. In p^] Constantin proved that the long-time average energy spectrum, E u LP (k) = 
limsup^^oo y Jq dt E^p(k, t) of the 2D Navier-Stokes solutions satisfies a bound of the form 



W LP (k)<C^k- 3 i^j , (2.16) 

where 7 = || Vu^H^oo. This upper bound is much larger than the Kraichnan-Batchelor spectrum 
( 2. ll| ) over most of the range ko < k < kd, but becomes comparable at the upper limit. We are 



going to assume here that something stronger is true of the 2D Navier-Stokes solutions, namely, 
for some T > 0: 

1 f T 

sup- / dt sup k 3 E£p(k,t) < 00 (2.17) 

Note that C u (t) = sup fc>fco k 3 E v LP {k , t) /rj 2 ^ is an instantaneous (worst) value of the Kraichnan- 
Batchelor constant and ( 2. 17] ) is a bound on its time-average. The most important aspect of 



this estimate — in contrast to what is so far proved, equation ( |2.16 ) — is its uniformity for small 
viscosity v > 0. Our main hypothesis ( |2.17| ) is equivalent to 

sup||u^|| L2 , 0T . B i,oo fT2 . n < 00 (2.18) 

Using estimates for singular integral operators, this may also be expressed equivalently in terms 
of the vorticity w" = Vxu", as 

sup II^II_l2( 0jT;B o, °°( T 2)) < 00 (2.19) 

This latter estimate could be stated in terms of the Littlewood-Paley enstrophy spectrum 
n u Lp (k,t) := k~ l \\uj v N {t)\\\ 2 for k G [2 N ,2 N+1 ), that 

1 f T 

sup- / dt sup kni P (k,t) < 00, (2.20) 

u>0 T Jq k>k 

which is entirely equivalent to the initial hypothesis (|2.17| ). 
We now state our third main result: 

Theorem 3 Let u u be the solution of the 2D incompressible Navier-Stokes equation for initial 
data vlq S B^ 00 ^ 2 ) and viscosity v > 0. Let uo v = Vxu 1 '. Assume that (2. It ) holds for these 
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solutions. Then, there exists a u G L 2 (0,T] B 2 '°° (Y 2 )) n Lip(0, T; i7~ L (T 2 )) *d 

solution of the 2D incompressible Euler equations in the velocity-pressure formulation, and for 

which, with u = Vxu, 



uj u ^uj weak-* in L 2 (0, T; B°'°°(T 2 )), (2.21) 



CO — > UJ 



strong in L 2 (0, T; V^" S ' 9 (T 2 )), (2.22) 



for some q > 2 and s > 1 — in the limit as v — > 0. Furthermore, u and u = K * l<j are a 
weafc solution of the 2D incompressible Euler equations in the vorticity-velocity formulation, in 
the sense that 

((d t + u-V)ip,u) =0 (&&s; 

for all tp G Cq°([0, T] x T 2 ) and i/ie expression ((<9< + u-V)^,w) is defined as the evaluation of 
a continuous linear functional on the element uj of the Banach space L 2 (0,T;B 2 '°°(T 2 )). 



This theorem essentially just states that the estimate ( |2.18| ) provides enough compactness to 
take limits along subsequences. Obviously, the hard problem is to prove that a bound such as 



( |2.18| ), as expected from 2D turbulence theory, really does hold. The theorem could be stated 
in a somewhat more general form, with the results on weak solutions in the velocity-pressure 
formulation remaining true for any p > 1 replacing p = 2, if a corresponding replacement is 



made in the estimate (2.1 8|) . Likewise, the results on weak solutions in the vorticity-velocity 



formulation will remain true for any p > 4/3 replacing p = 2. 

Our interest in this class of solutions is that they seem compatible with a finite rate of 
enstrophy dissipation in the inviscid limit. However, the very notion of "dissipative solution" 
must be reformulated. Local functions of the vorticity, of the form h(to(x,t)), do not need 
even to exist, since now the vorticity to is only a distribution and not necessarily a measurable 



function. Thus, a balance equation such as (L2) that we proved in Theorem 1 for DiPerna 



Majda solutions is not even well-defined for the class of solutions considered here. However, 
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the balance equations for the mollified vorticity in ( |2.2| ), namely, 

d t h(co E ) + V-[u e h(uj £ ) + h'(u e )tr e ] = -Z hj£ (u) (2.24) 

with Zf l£ (uj) = -/i"(w e )Vw £ '(T £ , are still perfectly well-defined. The term Zh e {uj) which appears 
as sink on the righthand side of (2.24) represents a flux of h to length-scales < e and it is expected 



to be asymptotically non-negative for small e. In fact, more should be true. A corresponding 
balance equation holds for the solutions of the 2D Navier-Stokes solutions uj u , in the form 

d t h{u v ) + V-[h(u u )u v - vVh(u v )} = -uh"{u v )\Vu) v \ 2 , (2.25) 

for any h S C 2 . Then we expect the following 

Conjecture 1 Let lo u be a sequence of solutions of the 2D Navier-Stokes equation obeying 
( 2.1$ ) and let u be the limiting 2D Euler solution, as provided by Theorem 3. Then, for this 



Euler solution 



Z k («) = limZ M H (2.26) 



exists in the sense of distributions for any h £ Ti.2- Furthermore, the same distribution is 
obtained by the limit of the viscous dissipation of the Navier-Stokes solutions: 

Z h (uj) = lim uh"{co u )\Vu; u \ 2 (2.27) 

v — >0 

for any h G Ti.2f]C 2 .In particular, for any such convex h, the distribution Zyjyu) is a nonneg- 
ative measure. Finally, there should exist a suitable such 2D Euler solution lo for which 

Z h (uj) > (2.28) 

with a strict inequality, for a convex h G H.2. 

The first limit statement in the conjecture may be put another way, which is perhaps more 
illuminating. Although the integral Ih(t) may itself be infinite for the Euler solutions in The- 
orem 3, it still makes sense to talk about a finite dissipation rate for it, defined as Dh(t) := 
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liminf £ _>o ~~3f{^)^ where I^(t) is the value of the integral for uj £ . The conjecture then states 



Dh(t) = J* T2 d 2 x Zf l {u)(yL,t) > 0. Note that for DiPerna-Majda solutions the first limit (2.26) 



has been demonstrated in Theorem 1 and it is easy to show for these solutions that the second 
limit ( gj^D also holds, using the same kind of argument as in Proposition 4 of Duchon-Robert 
|]]. Of course, for DiPerna-Majda solutions with u> G L p and p > 2 the distribution Z^iuS) = 
and thus the third statement is false. 

We believe that it is necessary to understand solutions of the type considered in Theorem 
3 in order to develop a rigorous mathematical theory of invariant measures for forced steady- 
states of 2D Navier-Stokes in the zero- viscosity limit. As proved in [jffl, Section 3.3.4, the mean 
enstrophy flux (Z e ) is a positive constant rj, independent of e, for length-scales e<^£f, the 
forcing scale, and e 3> zW 4 ^ 1 / 4 /77 1 / 4 , under the single assumption that the total mean energy 
E remains finite in the limit as v — > 0. (This requires adding an additional dissipation at 
low-wavenumbers to dispose of the "condensate" from the inverse energy cascade: see H). 
Thus, in the limit as v — > 0, we expect that the realizations of the ensemble shall be solutions 
of the (forced) 2D Euler equations with Z(u) > 0. If the statistical energy spectrum has 
the Batchelor-Kraichnan form in this limit, then individual realizations of the vorticity satisfy 
oj € £?2 _ '°°(T 2 ) a.s. It is not hard to prove this fact, using the methods of p|| (the wavelet 
characterization of Besov spaces and the Borel-Cantelli argument of Theorem 4). 

To see dissipation in the sense of our Conjecture 1 for the problem of free decay of 2D 
turbulence starting from random initial conditions, one must begin with initial data which 
is sufficiently rough. It is well-known that if one starts with loq € i?p'°°(T 2 ) for s > at 
time t = 0, then the exponent may (and generally will) deteriorate exponentially in time: for 
example, s(t) = e — Cll^oll^*^ if Uq e L°°(T 2 ) f| Bp' °°(T 2 ), but the exponent remains positive |||. 



Thus, there will be no dissipation at any finite time. On the other hand, the deterioration is 
consistent with the expectation from 2D turbulence theory that there will be an exponentially 



growing range of scales e with (Z £ (t)) « independent of e [^6|, ^7|. To see dissipation at 
finite (or zero) time, one must begin with initial data no more regular than loq G Bp'°°(T 2 ) a.s. 

13 



for p > 2. Such initial data could be prepared, for example, by taking an invariant measure for 
the driven problem and then turning off the force. The results of DiPerna and Lions [15] do 



not rule out dissipation in this instance, because they require u G L 1 (0, T; W 1,P (T 2 ) for some 
p > 1, whereas W l ^{T 2 ) = B^ P (T 2 ) C Sp'°°(T 2 ). If u(t) G Bp°°(T 2 ) only, then examples like 
that in section IV. 2 of |^] show that uniqueness of the Lagrangian trajectories breaks down 
and dissipation (in the sense of non-vanishing enstrophy flux) is possible. 

It is natural to expect that 2D Euler solutions which are dissipative in the proposed sense, 
i.e. Zh{uf) > for convex h, must be unique. Our Conjecture 1 states that "viscosity solutions" 
of 2D Euler equations are dissipative, so that these must then also be unique. Duchon and 
Robert have advanced the same idea for the 3D case. There is perhaps even more reason 
to believe so in 2D, because there is then an infinity of convex "entropies" h. For the problem 
of scalar conservation laws, such entropies play a crucial role in establishing uniqueness (e.g. 
see However, unlike the scalar case, it is not necessarily true even for smooth classical 

solutions of 2D Euler that the dynamics is L 1 -contractive. In fact, for two such solutions ui,u)2, 
Jj||k>i(t) — W2(i)lli = ~~ 2 f u =£J ni2*(ui — 112)0; ds where 1112 is the unit vector normal to the 
curve oj\ = u>2 from the region lu\ > u>2 to uj2 > and uj = oj\ = UJ2- It is precisely the nonlocal 
relation between u and u) which allows ui ^ 112 where u\ =0^2- So far, uniqueness of weak 
Euler solutions in 2D is established only for the solutions with u G L°°([0,T] xT 2 ) constructed 



by Yudovich [21] and for DiPerna-Majda solutions in L p ,p G (1, 00) if also u> G BMO [29]. It 
is not known in general whether the DiPerna-Majda solutions are unique, although for p > 2 
they are "dissipative Euler solutions" , in the sense that 

d t h{io) + V-[uh(u)] < (2.29) 

for all convex h G TC P . In fact, as noted above, DiPerna-Majda solutions for p > 2 are "renor- 
malized solutions" as considered by DiPerna-Lions and satisfy ( p.29| ) in the degenerate sense 
with equality. Yet their uniqueness is an open question. 
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3 Proofs 



2.1. Proof of Theorem 1 

We comment first on the validity of the weak vorticity- velocity equation for the DiPerna-Majda 
solutions. The condition p > 4/3 arises from the requirement that the nonlinear advection term 
uw G L 1 (T 2 ). Since u G W 1 ' P (T 2 ) C W' (T 2 ) for h = J - \ by Sobolev imbedding, one finds 
that > q, with g defined by - = 1 — -, when p > 4/3. Then uw G L 1 (T 2 ) follows by Holder 
inequality. The weak velocity-pressure form of the Euler equation is that 

J d 2 x J dt [d t cj)-M + V®</>:u®u]=0 (3.1) 

for any smooth, divergence-free test function </>(x, t). In particular, cf) = 'V ± ip satisfies these 
conditions for any smooth if}, where V -1- is the skew- gradient, = Sijdj with Eij the Levi- 
Civita tensor in 2D. (In fact, by Hodge theory, any divergence-free vector field cf) in 2D can be 
written in this way.) Substituting cf> = ^ ± ip into ( |3.l|) it is easy, using the L 1 property of uu; 
and u = — V^-u, to derive the vorticity- velocity equation by an approximation argument. 

The main condition of Theorem 1 on the index r can be similarly understood from the 
following lemma: 

Lemma 1 If lj £ L P (T 2 ) for p > 4/3 and u = K * uj, then for any h € TC r it holds that 
uh(oj) G L X (T 2 ) when r = |p — 1 for ^<p<2,r<p for p = 2, and r = p for p > 2. 

Remark: For convenience in the proof below, and in all later proofs, we employ an equivalent 
definition of the class of functions 

H p := {h\ h G C^R), \h'(u)\ < CM^ 1 for \to\ > R for some C,R > 0} (3.2) 

We will make the argument then assuming that R = so that the bound in ( |3.2| ) above holds 
globally. In fact, when R > it is easy to bound the contributions from the small-a; regions of 
integration over space and time by terms proportional to II/i'Uloo^^^, ll^'lli^f-^ij], assuming 
that the latter are finite. So we lose no generality and simplify the arguments by taking R = 0. 
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Proof of Lemma: We first note the definition K := V G where G is the Greens function 
of —A on T 2 . Then u E W l,p (Y 2 ) because ||u|| p < j|K||i ||cj|| p by Young's inequality and 
||Vu|| p < C||a;||p by the C alder on- Zygmund inequality. Hence, by the same Sobolev imbedding 
as before, u G L v> (T 2 ) for X = | — | when | < p < 2 and for any finite p' > 1 when p = 2, and 
for p' = oo when p > 2. Then, by definition of Tt r , 

\\uh(u))\\i < (const. )||u|cj| r ||i < (const. )||u|| p / ||u;||£ ? / (3-3) 

with 4 = 1 — A. When f < p < 2, then X = I — - and ro' = p for r = |p — 1. On the other 
hand, when p > 2, then </ = 1, and rg ; = p for r = p. Lastly, in the critical case p = 2, the 
only requirement is that q' > 1. Then rg' < p can be satisfied for any r < p by an appropriate 
choice of q' > 1. Thus, for the given definitions of r, 

||ufc(w)||i < (const.) (M^ 1 (3.4) 

because \\u\\ p > < C||u||^i,p < C"||ll>|| p . □ 

Proof of Theorem 1: We consider the filtered balance equation ( |2.2| ): 

d t h(u> £ ) + V-[u £ /i(w £ ) + fc,'(u; e )<r e ] = /i"(o%) Vuj £ -a £ . 

and, just as in Q, we show that every term on the lefthand side has a limit in the sense of 
distributions for e — > 0. We show first that h(uj £ ) — > /i(u>). In fact, by the mean-value theorem, 
/i(w £ ) — /i(u;) = h'{uj £ ){uj £ — to) for o) e (x, i) = A(x, i)u;(x, i) + (1 — A(x, i))u; e (x, t) with some 
< A(x, t) < 1. Then, in the notations of Lemma 1, we have 

\h(u £ (x,t)) — h(cv(x,t))\ < (const.) |a) e (x, t)\ r ~ l |w e (x, i) — w(x,t)| 

and thus by Holder inequality 

||/i(a; £ (i))-/j(u;(i))|| g , < (const.) ||u> £ (i) H^, 1 ||w e (t) - w(t)||r,' 

< (const.) \\uj(t)\\ r p - x \\u e {t)-uj(t) (3.5) 
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By the properties of the mollifier, lim^o ||^e(*) — = f° r a - e - t € [0, T], and thus 

lim £ ^o \\h(uJ £ (t)) ~~ h(co(t))\\ q i = 0. To complete the argument, we use the uniform bound 

\\h(u £ {t)) - h(u(t))\\ q , < (const.)||w||x,oo (0> T;LJ>(r«)) ( 3 - 6 ) 

to conclude by dominated convergence that lim £ ^o \\h(uJ e ) — h(uj)\\ Lq i^ Q T ] x j2) = 0, which implies 
convergence h(uj £ ) — > h(uj) in sense of distributions. 

We show next for the middle term that u £ h(u £ ) — > uh(u). In fact, with notations again as 
in Lemma 1, 

\\u £ (t)h(tO £ (t)) - U(t)h(uj(t))\\l < WMt) - "(t)\\p'\\HUe(t))\\ q > + ||u(t)||p,||/l(w e (t)) - h(to(t))\\ qf 

< (const.)K(i) - u(t)\\ p ,\\u;(t)\\ r p + \\n(t)\\ p/ \\h(co £ (t)) - h(co(t))\\ q , (3.7) 

Thus, we see that lim e ^o \\u £ h(uj £ (t)) — uh(oj(t))\\i = for a.e. t € [0, T\. In this case we have 
the uniform bound 

\\u £ h(u £ (t)) - Uh(w(t))||i < (conSt.)||w||2to 1 ( o,T;LP(Ta)) ( 3 ' 8 ) 

so that we can use Lebesgue's theorem again to infer lim e ^o || u £^(^e) — uh(u) IIl^Io.tJxT 2 ) = 0) 
which gives the result. 

Finally, for the third term we show that h'(u> £ )a £ -> as a distribution. We use the 
definition a £ = (uto) £ — u £ u> £ and the Holder inequality 

\\h'(u e (t))<r e (t)\\i < ll^(w e (t))|| P /(r-i)||(u(t)w(t)) e -u £ (t)u £ (t)\\ p/(p _ r+1) (3.9) 
along with \\h' (uj £ (t))\\ p ^ r _ 1 ^ < (const.) ||w e (t) 1 1 ^ 1 and the triangle inequality 

||(u(t)w(t)) e -u £ (t)u £ (t)\\ p/(j) _ r+1) < \\(u{t)u(t)) £ - u(t)u{t)\\ p/{p _ r+1) 

+ ||u(t)|| l/ ||w(t)-w e (t)|| I , + ||u(t)-u e (t)|| p /||w e (t)|| p (3.10) 

to infer that lim e ^o \\{ u (t)uj(t)) £ — u £ (t)u £ (t) ||p/(p— r-+i) = for a.e. t <G [0, T]. Note that we 
have used \\u(t)u(t)\\ p /^ p _ r+ i^ < \\u(t)\\ p /( p _ r )\\ui(t)\\ p and (p — r)/p < p'. Again a uniform 
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bound on \\h' (io £ (t))cr £ (t)\\i like that in ( |3.8| ) completes the argument. Gathering these results, 
we see that the entire lefthand side of ( |2.2| ) approaches dth(u) + V-[u/i(w)] in the sense of 
distributions as e — > 0. Obviously this limit is independent of the mollifier ip and the righthand 
side — Z]- he {uj) = h" (uj £ )'Vuj £ -cr £ has the same limit. This gives the first half of Theorem 1. 

The second half of the theorem for the particular choice h(uj) = ^|w| 2 follows by the same 
argument as in |j. In this proof, the balance ( |2.2| ) is replaced by 

d t (~UU £ ) + V-[(±UM e )u] = ~Z e (u) (3.11) 

where an easy calculation gives 

Z £ {uj) := io;V-[(a;u) e ] - ~cj(u-V)u £ . (3.12) 

An argument exactly like the previous one shows that, when p > 2, the distributional limit 
lim e _+o Z £ (lo) exists and equals —dt(^\uj\ 2 ) — V-[(||u;| 2 )u] = Z(u>). In addition, a simple calcu- 
lation using the incompressibility of the velocity field shows that the expression appearing in 
( |2.7| ) in Theorem 1 can be written 

f dH V(^)-AHAH 2 = V-[u(^ 2 ) e - (uw 2 ) £ ] + 4Z £ (lo). (3.13) 

As before, it is easy to show for p > 2 that u(u; 2 ) e — (uuj 2 ) £ — > as a distribution when e — > 0. 
Hence, it follows that the limits of \ f d 2 £ V^ £ (^)-A^u| A^w| 2 and Z £ (uj) are also the same. 
That proves the second half of Theorem 1. □ 
2.2. Proof of Theorem 2 

A result on global conservation corresponding to the local result in Theorem 2 was already 
proved in || but with an additional smoothness assumption that u G L p (0, T; i?p'°°(T 2 )). Here 
we show that conservation holds without any such a smoothness assumption. Let T £ (f,g) := 
(fd)e ~ fs9s where f £ = <p £ * f ■ Then, we make use of the following key estimate: 

Lemma 2 Let u G L P (T 2 ) and u G W^'^T 2 ) for p>2, and let V-u = 0. Then 

||V-T £ (u,u;)|| L p/2 < C||u||wi,p|M|.&p (3-U) 
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with a constant C independent of e. 
Proof: Note that 

V-T e (u,a>) = V-[(uw) E - ULU £ ] + (u - u £ )-Vlo £ . (3.15) 

The first term is handled in exactly the same manner as in Lemma II. 1 of [15]. However, it is 
easy to see that 

|| u — u £ ||lp < e||Vu||iP < e||u||^i,p (3.16) 

and 

||Vu; e ||i,p < e^ 1 1| V^||x,i ||u;||iP. (3-17) 
These control the second term. □ 

Corollary 1 Under the same hypotheses, let r £ := — V-r e (u, uj). Then lim e ^o r £ = strong in 
Wl 2 (T 2 ) forp>2. 

Proof: Since lim e ^o r e = for smooth u, cj, one can obtain the result for all to S L p (T 2 ),u € 
W 1,P (T 2 ) by an approximation argument using the estimate in Proposition 1. □ 

If u is related to u by the Biot-Savart formula, u = K * u, then r e (u, u) = cr £ in the earlier 
notation. In particular, we see that 

8 t uj £ + (iWK = r £ (3.18) 

for a weak Euler solution. 

Proof of Theorem 2: Using (|3.1Sj) we get 



dth(u e ) + (u £ -V)h(u £ ) = h'(uj £ )r £ . (3.19) 
It was proved in Theorem 1 that 

d t h(uj £ ) + (u £ -V)/i(tj e ) — ► d t h(io) + (u-V)/i(w) (3.20) 
in the sense of distributions for all such h. Furthermore, for any h € C 1 with h' G L°°, 

\\ti{u s )r s \\ L i^0. (3.21) 
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Having proved that (2.9) holds for h with h' 6 L°° we then extend it to the general h in the 



theorem statement by an approximation argument, as in Corollaries II. 1-2 in [15]. □ 
Remark: The smoothness assumed in the earlier proof of [J7J is not necessary to obtain con- 
servation, but only to provide an estimate of the rate of the vanishing of the flux. With the 
assumption that to € L p (0,T; Bp ,(X) (T 2 )) the bounds above can be improved as follows. Gen- 
eral estimates in Besov spaces give ||Vu; e (t)||p < Ce s_1 ||a;(i)|| B s,oo and sup^| <£ ||A^u;(i)||p < 
e s ||u;(i)|| B s,°°. See fl2|| , or Appendix C of j7]]. Just as in Q this gives 



\\Z h 4u(t))\\i < (constVl^H^IKt)!!^ < (const.yi^i)!^ (3.22) 

Because r + 1 < p, integrating over t € [0, T] gives 

ll^, E (^)llxi([0,T]xT2) < (const.)£ 2s ||a;||^+ ( 1 0>T . B »,<x. (T2)) . (3.23) 

Thus, lim £ ^o %h e(^) = as before, but with an estimate of the rate. The bound 0(e 2s ) is in 
agreement with the estimate given by the heuristic argument in the Introduction. 

2.3. Proof of Theorem 3 

Theorem 3 is a consequence of the following technical lemma: 

Proposition 1 Consider a sequence {uj £ \e > 0} and u £ = K * co £ given by the Biot-Savart 
formula, with the following properties: 

sup||u; e || Lr(0iT . B o,oc (T2)) < oo (3.24) 



for r,p € [2, oo], and 



SUp||u £ ||Lip(0,T;H-i(T2)) < OO (3.25) 

£>0 



for some L > 3. Then, there exist uj and u = K * uj with 

uj G L r (0,T; Bp ,oc (T 2 )) (3.26) 

and 

u € Lip(0,r;^r L (T 2 )). (3.27) 
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and, furthermore, there exists a subsequence ofu> £ ,u £ along which 



oj £ ^u weak-* in L r (0, T; B® ,00 (T 2 )), 



U) — ► OJ 



strong in L r (0,T;W- s > q (T 2 )), 



(3.28) 
(3.29) 



for some q > p and s > 2 ( i — ^ ) and for t = min{r, q} > 2, 



u £ ^u strong in L*([0, T] x T 2 ). 



Proof: The first statement ( [3.28 ) on weak-* convergence of u> £ to uj € L r (0, T; 5p'°°(T 2 )) is a 
simple consequence of the Banach-Alaoglu theorem. 



We derive the second statement from the Aubin-Lions compactness criterion (see [ 30 1 



Theorem 5.1 or [31], Theorem III. 2.1). Note first that there is the continuous embedding 
Bp'°°(T 2 ) C Bq S '' q (T 2 ) = W- s '' g (T 2 ) for each q > p and s' > 2 (± - ±) (see @, Theorem 



2.7.1 and Prop. 2.3. 2/2). Therefore, from (|3.24|) , 



SU P ll^ 6 IIl , '(0,T;VF- s '''?(T 2 )) ^ 00 



£>0 



(3.31) 



On the other hand, from ( |3.25| ), 



sup 

£>0 



dw £ 



dt 



L°°(0,T;H-( L + 1 )(T 2 )) 



(3.32) 



Furthermore, for s > s 1 and L + 1 > s there are continuous embeddings 



W - S ', q( j2^ c ^--.g^jj c H- {L+1) (T 2 ), 



(3.33) 



and the first embedding is compact by the Rellich-Kondrachov theorem (see |p2l , Chapter 12). 
Hence, we conclude that {ui £ \e > 0} is compact in -L r (0, T; W~ s,q (T 2 )) and contains a strongly 
convergent subsequence. 

To obtain the third result we remark that one may choose < s < 1 and that the mapping 
uj \— >u = K*iijis continuous from ^"^(T 2 ) into W 1 ~ s ' q (T 2 ), because of the continuity of the 
singular integral operator T(cj) = (VK) * uj from W~ s,q (T 2 ) into itself (for example, see |33| , 
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Theorem 3.2.1) and the bound HuH^i-s,?^) < (const.) [||u||w-s,«(i'a) + || Vu|| w - S , 9 ( T 2)] (see 
[12], Theorem 2.3.8). Of course, convergence of u £ — ► u strong in L r (0,T; M /1_S,,? (T 2 )) implies 
at once convergence strong in L*([0,T] x T 2 ). □ 

Proof of Theorem 3: The proof is very straightforward and quite similar to that of DiPerna and 
Majda in |ll| for ujq € L p with p > 2 (the easier case than 1 < p < 2). In fact, the Lipschitz 
estimate in time 

sup||u I/ || Ijip(0iT . jff - £(T 2 )) < oo (3.34) 

holds for the 2D Leray solutions with initial energy finite, Eq := ^||uo||| 2 < °o, which is part 
of our assumption. See Section 2A and Appendix A of [11], for example. But, in that case, 
the Proposition 1 applies, with r = p = 2. The limiting velocity u is easily seen to be a weak 
solution of the 2D Euler equation in the velocity-pressure formulation, because of the third 
result ( |3.30| ), the strong L* convergence w v — > u with t > 2. Obviously, more general versions of 
Theorem 3 for any r,p£ [2, oo] could be proved, with (3.24) replacing ( p.!9| ) in the hypothesis. 

The statements on the weak solutions in the vorticity-velocity formulation follow from ar- 
guments very similar to those earlier in Lemma 1, but now using the density of C°°(T 2 ) in 
Bq' 1 ^ 2 ) Theorem 2.3.3). So, we just verify the required regularity of u. Because of the 

hypothesis on u), u € L r (0,T; Bp ,oc (T 2 )). Indeed, using the C alder on- Zygmund inequality it is 
easy to show that Vu G L r (0, T; Bp ,co (T 2 )), and this is equivalent to the first statement ([12], 
Theorem 2.3.8). Then, for any p" < p' with p' defined by ^ = - — |, one has the continuous 
embedding Bp°°(T 2 ) C B° pl }(T 2 ). In fact, ^'^(T 2 ) C B £ p ,?°(T 2 ) for ^ := £ + f for any small 
e > ((HI, Theorem 2.7.1), but then Bpf°(T 2 ) C B^(T 2 ) by an elementary imbedding (p|, 
Prop. 2.3.2/2). Now, precisely for p > |, one has \ < | with ^ + - = 1. Thus, it is possible 
to choose p" > q but still p" < p'. Also, with t defined by - + i = 1, r > t for r > 2. In that 
case, u G L*(0, T; Sg'^T 2 )), whereas L r (0, T; 5p'°°(T 2 )) = [L*(0, T; ^^(T 2 ))]*, the Banach 
dual (||], Theorem 2.11.2). To conclude the proof, we just note that, if ip G C£°([0,T] x T 2 ), 
then also (d t + u-V)V> G L*(0, T; 5° ,:L (T 2 )) (see ||], Lemma 3.3.1). □ 
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